DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs.

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Activation of Notch and Myc Signaling via B-cell-Restricted Depletion of Dnmt3a Generates a Consistent Murine Model of Chronic Lymphocytic Leukemia.

Chronic lymphocytic leukemia (CLL) is characterized by disordered DNA methylation, suggesting these epigenetic changes might play a critical role in disease onset and progression. The methyltransferase DNMT3A is a key regulator of DNA methylation. Although DNMT3A somatic mutations in CLL are rare, we found that low DNMT3A expression is associated with more aggressive disease. A conditional knockout mouse model showed that homozygous depletion of Dnmt3a from B cells results in the development of CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yields a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haploinsufficient tumor suppressor. B1a cells were confirmed as the cell of origin of disease in this model, and Dnmt3a depletion resulted in focal hypomethylation and activation of Notch and Myc signaling. Amplification of chromosome 15 containing the Myc gene was detected in all CLL mice tested, and infiltration of high-Myc-expressing CLL cells in the spleen was observed. Notably, hyperactivation of Notch and Myc signaling was exclusively observed in the Dnmt3a CLL mice, but not in three other CLL mouse models tested (Sf3b1-Atm, Ikzf3, and MDR), and Dnmt3a-depleted CLL were sensitive to pharmacologic inhibition of Notch signaling in vitro and in vivo. Consistent with these findings, human CLL samples with lower DNMT3A expression were more sensitive to Notch inhibition than those with higher DNMT3A expression. Altogether, these results suggest that Dnmt3a depletion induces CLL that is highly dependent on activation of Notch and Myc signaling. SIGNIFICANCE: Loss of DNMT3A expression is a driving event in CLL and is associated with aggressive disease, activation of Notch and Myc signaling, and enhanced sensitivity to Notch inhibition.

Smart-RRBS for single-cell methylome and transcriptome analysis.

The integration of DNA methylation and transcriptional state within single cells is of broad interest. Several single-cell dual- and multi-omics approaches have been reported that enable further investigation into cellular heterogeneity, including the discovery and in-depth study of rare cell populations. Such analyses will continue to provide important mechanistic insights into the regulatory consequences of epigenetic modifications. We recently reported a new method for profiling the DNA methylome and transcriptome from the same single cells in a cancer research study. Here, we present details of the protocol and provide guidance on its utility. Our Smart-RRBS (reduced representation bisulfite sequencing) protocol combines Smart-seq2 and RRBS and entails physically separating mRNA from the genomic DNA. It generates paired epigenetic promoter and RNA-expression measurements for ~24% of protein-coding genes in a typical single cell. It also works for micro-dissected tissue samples comprising hundreds of cells. The protocol, excluding flow sorting of cells and sequencing, takes ~3 d to process up to 192 samples manually. It requires basic molecular biology expertise and laboratory equipment, including a PCR workstation with UV sterilization, a DNA fluorometer and a microfluidic electrophoresis system.

Discovery of Candidate DNA Methylation Cancer Driver Genes.

Epigenetic alterations, such as promoter hypermethylation, may drive cancer through tumor suppressor gene inactivation. However, we have limited ability to differentiate driver DNA methylation (DNAme) changes from passenger events. We developed DNAme driver inference-MethSig-accounting for the varying stochastic hypermethylation rate across the genome and between samples. We applied MethSig to bisulfite sequencing data of chronic lymphocytic leukemia (CLL), multiple myeloma, ductal carcinoma in situ, glioblastoma, and to methylation array data across 18 tumor types in TCGA. MethSig resulted in well-calibrated quantile-quantile plots and reproducible inference of likely DNAme drivers with increased sensitivity/specificity compared with benchmarked methods. CRISPR/Cas9 knockout of selected candidate CLL DNAme drivers provided a fitness advantage with and without therapeutic intervention. Notably, DNAme driver risk score was closely associated with adverse outcome in independent CLL cohorts. Collectively, MethSig represents a novel inference framework for DNAme driver discovery to chart the role of aberrant DNAme in cancer. SIGNIFICANCE: MethSig provides a novel statistical framework for the analysis of DNA methylation changes in cancer, to specifically identify candidate DNA methylation driver genes of cancer progression and relapse, empowering the discovery of epigenetic mechanisms that enhance cancer cell fitness.This article is highlighted in the In This Issue feature, p. 2113.

Preneoplastic Alterations Define CLL DNA Methylome and Persist through Disease Progression and Therapy.

Most human cancers converge to a deregulated methylome with reduced global levels and elevated methylation at select CpG islands. To investigate the emergence and dynamics of the cancer methylome, we characterized genome-wide DNA methylation in pre-neoplastic monoclonal B cell lymphocytosis (MBL) and chronic lymphocytic leukemia (CLL), including serial samples collected across disease course. We detected the aberrant tumor-associated methylation landscape at CLL diagnosis and found no significantly differentially methylated regions in the high-count MBL-to-CLL transition. Patient methylomes showed remarkable stability with natural disease and post-therapy progression. Single CLL cells were consistently aberrantly methylated, indicating a homogeneous transition to the altered epigenetic state, and a distinct expression profile together with MBL cells compared to normal B cells. Our longitudinal analysis reveals the cancer methylome to emerge early, which may provide a platform for subsequent genetically-driven growth dynamics and together with its persistent presence suggests a central role in the normal-to-cancer transition.

Distinct evolutionary paths in chronic lymphocytic leukemia during resistance to the graft-versus-leukemia effect.

Leukemic relapse remains a major barrier to successful allogeneic hematopoietic stem cell transplantation (allo-HSCT) for aggressive hematologic malignancies. The basis for relapse of advanced lymphoid malignancies remains incompletely understood and may involve escape from the graft-versus-leukemia (GvL) effect. We hypothesized that for patients with chronic lymphocytic leukemia (CLL) treated with allo-HSCT, leukemic cell-intrinsic features influence transplant outcomes by directing the evolutionary trajectories of CLL cells. Integrated genetic, transcriptomic, and epigenetic analyses of CLL cells from 10 patients revealed that the clinical kinetics of post-HSCT relapse are shaped by distinct molecular dynamics. Early relapses after allo-HSCT exhibited notable genetic stability; single CLL cell transcriptional analysis demonstrated a cellular heterogeneity that was static over time. In contrast, CLL cells relapsing late after allo-HSCT displayed notable genetic evolution and evidence of neoantigen depletion, consistent with marked single-cell transcriptional shifts that were unique to each patient. We observed a greater rate of epigenetic change for late relapses not seen in early relapses or relapses after chemotherapy alone, suggesting that the selection pressures of the GvL bottleneck are unlike those imposed by chemotherapy. No selective advantage for human leukocyte antigen (HLA) loss was observed, even when present in pretransplant subpopulations. Gain of stem cell modules was a common signature associated with leukemia relapse regardless of posttransplant relapse kinetics. These data elucidate the biological pathways that underlie GvL resistance and posttransplant relapse.

TETs compete with DNMT3 activity in pluripotent cells at thousands of methylated somatic enhancers.

Mammalian cells stably maintain high levels of DNA methylation despite expressing both positive (DNMT3A/B) and negative (TET1-3) regulators. Here, we analyzed the independent and combined effects of these regulators on the DNA methylation landscape using a panel of knockout human embryonic stem cell (ESC) lines. The greatest impact on global methylation levels was observed in DNMT3-deficient cells, including reproducible focal demethylation at thousands of normally methylated loci. Demethylation depends on TET expression and occurs only when both DNMT3s are absent. Dynamic loci are enriched for hydroxymethylcytosine and overlap with subsets of putative somatic enhancers that are methylated in ESCs and can be activated upon differentiation. We observe similar dynamics in mouse ESCs that were less frequent in epiblast stem cells (EpiSCs) and scarce in somatic tissues, suggesting a conserved pluripotency-linked mechanism. Taken together, our data reveal tightly regulated competition between DNMT3s and TETs at thousands of somatic regulatory sequences within pluripotent cells.

Small Molecule Targets TMED9 and Promotes Lysosomal Degradation to Reverse Proteinopathy.

Intracellular accumulation of misfolded proteins causes toxic proteinopathies, diseases without targeted therapies. Mucin 1 kidney disease (MKD) results from a frameshift mutation in the MUC1 gene (MUC1-fs). Here, we show that MKD is a toxic proteinopathy. Intracellular MUC1-fs accumulation activated the ATF6 unfolded protein response (UPR) branch. We identified BRD4780, a small molecule that clears MUC1-fs from patient cells, from kidneys of knockin mice and from patient kidney organoids. MUC1-fs is trapped in TMED9 cargo receptor-containing vesicles of the early secretory pathway. BRD4780 binds TMED9, releases MUC1-fs, and re-routes it for lysosomal degradation, an effect phenocopied by TMED9 deletion. Our findings reveal BRD4780 as a promising lead for the treatment of MKD and other toxic proteinopathies. Generally, we elucidate a novel mechanism for the entrapment of misfolded proteins by cargo receptors and a strategy for their release and anterograde trafficking to the lysosome.

Epigenetic evolution and lineage histories of chronic lymphocytic leukaemia.

Genetic and epigenetic intra-tumoral heterogeneity cooperate to shape the evolutionary course of cancer1. Chronic lymphocytic leukaemia (CLL) is a highly informative model for cancer evolution as it undergoes substantial genetic diversification and evolution after therapy2,3. The CLL epigenome is also an important disease-defining feature4,5, and growing populations of cells in CLL diversify by stochastic changes in DNA methylation known as epimutations6. However, previous studies using bulk sequencing methods to analyse the patterns of DNA methylation were unable to determine whether epimutations affect CLL populations homogeneously. Here, to measure the epimutation rate at single-cell resolution, we applied multiplexed single-cell reduced-representation bisulfite sequencing to B cells from healthy donors and patients with CLL. We observed that the common clonal origin of CLL results in a consistently increased epimutation rate, with low variability in the cell-to-cell epimutation rate. By contrast, variable epimutation rates across healthy B cells reflect diverse evolutionary ages across the trajectory of B cell differentiation, consistent with epimutations serving as a molecular clock. Heritable epimutation information allowed us to reconstruct lineages at high-resolution with single-cell data, and to apply this directly to patient samples. The CLL lineage tree shape revealed earlier branching and longer branch lengths than in normal B cells, reflecting rapid drift after the initial malignant transformation and a greater proliferative history. Integration of single-cell bisulfite sequencing analysis with single-cell transcriptomes and genotyping confirmed that genetic subclones mapped to distinct clades, as inferred solely on the basis of epimutation information. Finally, to examine potential lineage biases during therapy, we profiled serial samples during ibrutinib-associated lymphocytosis, and identified clades of cells that were preferentially expelled from the lymph node after treatment, marked by distinct transcriptional profiles. The single-cell integration of genetic, epigenetic and transcriptional information thus charts the lineage history of CLL and its evolution with therapy.

Corrupted coordination of epigenetic modifications leads to diverging chromatin states and transcriptional heterogeneity in CLL.

Cancer evolution is fueled by epigenetic as well as genetic diversity. In chronic lymphocytic leukemia (CLL), intra-tumoral DNA methylation (DNAme) heterogeneity empowers evolution. Here, to comprehensively study the epigenetic dimension of cancer evolution, we integrate DNAme analysis with histone modification mapping and single cell analyses of RNA expression and DNAme in 22 primary CLL and 13 healthy donor B lymphocyte samples. Our data reveal corrupted coherence across different layers of the CLL epigenome. This manifests in decreased mutual information across epigenetic modifications and gene expression attributed to cell-to-cell heterogeneity. Disrupted epigenetic-transcriptional coordination in CLL is also reflected in the dysregulation of the transcriptional output as a function of the combinatorial chromatin states, including incomplete Polycomb-mediated gene silencing. Notably, we observe unexpected co-mapping of typically mutually exclusive activating and repressing histone modifications, suggestive of intra-tumoral epigenetic diversity. Thus, CLL epigenetic diversification leads to decreased coordination across layers of epigenetic information, likely reflecting an admixture of cells with diverging cellular identities.

Targets and genomic constraints of ectopic Dnmt3b expression.

DNA methylation plays an essential role in mammalian genomes and expression of the responsible enzymes is tightly controlled. Deregulation of the de novo DNA methyltransferase DNMT3B is frequently observed across cancer types, yet little is known about its ectopic genomic targets. Here, we used an inducible transgenic mouse model to delineate rules for abnormal DNMT3B targeting, as well as the constraints of its activity across different cell types. Our results explain the preferential susceptibility of certain CpG islands to aberrant methylation and point to transcriptional state and the associated chromatin landscape as the strongest predictors. Although DNA methylation and H3K27me3 are usually non-overlapping at CpG islands, H3K27me3 can transiently co-occur with DNMT3B-induced DNA methylation. Our genome-wide data combined with ultra-deep locus-specific bisulfite sequencing suggest a distributive activity of ectopically expressed Dnmt3b that leads to discordant CpG island hypermethylation and provides new insights for interpreting the cancer methylome.

Cancer-Germline Antigen Expression Discriminates Clinical Outcome to CTLA-4 Blockade.

CTLA-4 immune checkpoint blockade is clinically effective in a subset of patients with metastatic melanoma. We identify a subcluster of MAGE-A cancer-germline antigens, located within a narrow 75 kb region of chromosome Xq28, that predicts resistance uniquely to blockade of CTLA-4, but not PD-1. We validate this gene expression signature in an independent anti-CTLA-4-treated cohort and show its specificity to the CTLA-4 pathway with two independent anti-PD-1-treated cohorts. Autophagy, a process critical for optimal anti-cancer immunity, has previously been shown to be suppressed by the MAGE-TRIM28 ubiquitin ligase in vitro. We now show that the expression of the key autophagosome component LC3B and other activators of autophagy are negatively associated with MAGE-A protein levels in human melanomas, including samples from patients with resistance to CTLA-4 blockade. Our findings implicate autophagy suppression in resistance to CTLA-4 blockade in melanoma, suggesting exploitation of autophagy induction for potential therapeutic synergy with CTLA-4 inhibitors.

Global delay in nascent strand DNA methylation.

Cytosine methylation is widespread among organisms and essential for mammalian development. In line with early postulations of an epigenetic role in gene regulation, symmetric CpG methylation can be mitotically propagated over many generations with extraordinarily high fidelity. Here, we combine BrdU labeling and immunoprecipitation with genome-wide bisulfite sequencing to explore the inheritance of cytosine methylation onto newly replicated DNA in human cells. Globally, we observe a pronounced lag between the copying of genetic and epigenetic information in embryonic stem cells that is reconsolidated within hours to accomplish faithful mitotic transmission. Populations of arrested cells show a global reduction of lag-induced intermediate CpG methylation when compared to proliferating cells, whereas sites of transcription factor engagement appear cell-cycle invariant. Alternatively, the cancer cell line HCT116 preserves global epigenetic heterogeneity independently of cell-cycle arrest. Taken together, our data suggest that heterogeneous methylation largely reflects asynchronous proliferation, but is intrinsic to actively engaged cis-regulatory elements and cancer.

Genetic determinants and epigenetic effects of pioneer-factor occupancy.

Transcription factors (TFs) direct developmental transitions by binding to target DNA sequences, influencing gene expression and establishing complex gene-regultory networks. To systematically determine the molecular components that enable or constrain TF activity, we investigated the genomic occupancy of FOXA2, GATA4 and OCT4 in several cell types. Despite their classification as pioneer factors, all three TFs exhibit cell-type-specific binding, even when supraphysiologically and ectopically expressed. However, FOXA2 and GATA4 can be distinguished by low enrichment at loci that are highly occupied by these factors in alternative cell types. We find that expression of additional cofactors increases enrichment at a subset of these sites. Finally, FOXA2 occupancy and changes to DNA accessibility can occur in G1-arrested cells, but subsequent loss of DNA methylation requires DNA replication.

Epigenetic restriction of extraembryonic lineages mirrors the somatic transition to cancer.

In mammals, the canonical somatic DNA methylation landscape is established upon specification of the embryo proper and subsequently disrupted within many cancer types. However, the underlying mechanisms that direct this genome-scale transformation remain elusive, with no clear model for its systematic acquisition or potential developmental utility. Here, we analysed global remethylation from the mouse preimplantation embryo into the early epiblast and extraembryonic ectoderm. We show that these two states acquire highly divergent genomic distributions with substantial disruption of bimodal, CpG density-dependent methylation in the placental progenitor. The extraembryonic epigenome includes specific de novo methylation at hundreds of embryonically protected CpG island promoters, particularly those that are associated with key developmental regulators and are orthologously methylated across most human cancer types. Our data suggest that the evolutionary innovation of extraembryonic tissues may have required co-option of DNA methylation-based suppression as an alternative to regulation by Polycomb-group proteins, which coordinate embryonic germ-layer formation in response to extraembryonic cues. Moreover, we establish that this decision is made deterministically, downstream of promiscuously used-and frequently oncogenic-signalling pathways, via a novel combination of epigenetic cofactors. Methylation of developmental gene promoters during tumorigenesis may therefore reflect the misappropriation of an innate trajectory and the spontaneous reacquisition of a latent, developmentally encoded epigenetic landscape.

Prolonged Mek1/2 suppression impairs the developmental potential of embryonic stem cells.

Concomitant activation of the Wnt pathway and suppression of Mapk signalling by two small molecule inhibitors (2i) in the presence of leukaemia inhibitory factor (LIF) (hereafter termed 2i/L) induces a naive state in mouse embryonic stem (ES) cells that resembles the inner cell mass (ICM) of the pre-implantation embryo. Since the ICM exists only transiently in vivo, it remains unclear how sustained propagation of naive ES cells in vitro affects their stability and functionality. Here we show that prolonged culture of male mouse ES cells in 2i/L results in irreversible epigenetic and genomic changes that impair their developmental potential. Furthermore, we find that female ES cells cultured in conventional serum plus LIF medium phenocopy male ES cells cultured in 2i/L. Mechanistically, we demonstrate that the inhibition of Mek1/2 is predominantly responsible for these effects, in part through the downregulation of DNA methyltransferases and their cofactors. Finally, we show that replacement of the Mek1/2 inhibitor with a Src inhibitor preserves the epigenetic and genomic integrity as well as the developmental potential of ES cells. Taken together, our data suggest that, although short-term suppression of Mek1/2 in ES cells helps to maintain an ICM-like epigenetic state, prolonged suppression results in irreversible changes that compromise their developmental potential.

Development and Validation of a Mass Spectrometry-Based Assay for the Molecular Diagnosis of Mucin-1 Kidney Disease.

Mucin-1 kidney disease, previously described as medullary cystic kidney disease type 1 (MCKD1, OMIM 174000), is an autosomal dominant tubulointerstitial kidney disease recently shown to be caused by a single-base insertion within the variable number tandem repeat region of the MUC1 gene. Because of variable age of disease onset and often subtle signs and symptoms, clinical diagnosis of mucin-1 kidney disease and differentiation from other forms of hereditary kidney disease have been difficult. The causal insertion resides in a variable number tandem repeat region with high GC content, which has made detection by standard next-generation sequencing impossible to date. The inherently difficult nature of this mutation required an alternative method for routine detection and clinical diagnosis of the disease. We therefore developed and validated a mass spectrometry-based probe extension assay with a series of internal controls to detect the insertion event using 24 previously characterized positive samples from patients with mucin-1 kidney disease and 24 control samples known to be wild type for the variant. Validation results indicate an accurate and reliable test for clinically establishing the molecular diagnosis of mucin-1 kidney disease with 100% sensitivity and specificity across 275 tests called.

Locally disordered methylation forms the basis of intratumor methylome variation in chronic lymphocytic leukemia.

Intratumoral heterogeneity plays a critical role in tumor evolution. To define the contribution of DNA methylation to heterogeneity within tumors, we performed genome-scale bisulfite sequencing of 104 primary chronic lymphocytic leukemias (CLLs). Compared with 26 normal B cell samples, CLLs consistently displayed higher intrasample variability of DNA methylation patterns across the genome, which appears to arise from stochastically disordered methylation in malignant cells. Transcriptome analysis of bulk and single CLL cells revealed that methylation disorder was linked to low-level expression. Disordered methylation was further associated with adverse clinical outcome. We therefore propose that disordered methylation plays a similar role to that of genetic instability, enhancing the ability of cancer cells to search for superior evolutionary trajectories.

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells.

Recent molecular studies have shown that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels and phenotypic output, with important functional consequences. Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs or proteins simultaneously, because genomic profiling methods could not be applied to single cells until very recently. Here we use single-cell RNA sequencing to investigate heterogeneity in the response of mouse bone-marrow-derived dendritic cells (BMDCs) to lipopolysaccharide. We find extensive, and previously unobserved, bimodal variation in messenger RNA abundance and splicing patterns, which we validate by RNA-fluorescence in situ hybridization for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, we identify a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, we show that variability in this module may be propagated through an interferon feedback circuit, involving the transcriptional regulators Stat2 and Irf7. Our study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.

Mutations causing medullary cystic kidney disease type 1 lie in a large VNTR in MUC1 missed by massively parallel sequencing.

Although genetic lesions responsible for some mendelian disorders can be rapidly discovered through massively parallel sequencing of whole genomes or exomes, not all diseases readily yield to such efforts. We describe the illustrative case of the simple mendelian disorder medullary cystic kidney disease type 1 (MCKD1), mapped more than a decade ago to a 2-Mb region on chromosome 1. Ultimately, only by cloning, capillary sequencing and de novo assembly did we find that each of six families with MCKD1 harbors an equivalent but apparently independently arising mutation in sequence markedly under-represented in massively parallel sequencing data: the insertion of a single cytosine in one copy (but a different copy in each family) of the repeat unit comprising the extremely long (∼1.5-5 kb), GC-rich (>80%) coding variable-number tandem repeat (VNTR) sequence in the MUC1 gene encoding mucin 1. These results provide a cautionary tale about the challenges in identifying the genes responsible for mendelian, let alone more complex, disorders through massively parallel sequencing.